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Abstract. Let f(n) denote the number of odd entries in the n th 
row of Pascal's binomial triangle. We study "average dispersion" and "typical 
dispersion" of f(n) - the latter involves computing a generalized Lyapunov 
exponent - and then turn to numerical analysis of higher dimensional examples. 

Let D and Di denote m x m nonnegative matrices. Let Zq, Z\, . . ., z k-i 
denote a sequence of independent random coin tosses (heads=l and tails=0 with 
equal probability). The Lyapunov exponent corresponding to random products of D Q 
and D% is 

1 



A 



lim — In 



D Z0 D Zl ■ ■ ■ D Zk2 D Zk _ 1 



almost surely. We computed A in an earlier paper pQ for a number of examples; our 
purpose was to quantify the "typical growth" of certain number-theoretic functions. 
To assess the corresponding "typical dispersion", a second-order extension of A is 
needed. 

The generalized Lyapunov exponent or moment Lyapunov exponent is 



L(t) 



lim — In ( E 

fc^oo k 



D z ,D zl ...D Zh _ 2 D Zh _ x 



for real t. Clearly 1/(0) = 0, that is, e L ^ = 1. Also, e L ^ is the largest eigenvalue in 
modulus of | (Do + D±), and e L ^ is the largest eigenvalue of \(Dq ® D + D\® Di), 
where <8> is the direct or Kronecker product of matrices. The latter result is called 
the replica trick |2J, which can be applied for arbitrary integer t > 3 as well. 

For simplicity, write z = z z 1 . . . z k _ 2 Zk-i and D z = D Z0 D Zl ■ ■ ■ D Zk 2 D Zk v Dif- 
ferentiating L(t) with respect to t, we obtain 



L'(0)= lim-E(ln(||£>,||)) = A > 



L"(0) 



lim — 

fc^oo k 



E (in i 



\D, 



E(ln(||D 2 



lim - Var (ln(||D 2 

k— >oo k 



a 



This definition of a 2 , unfortunately, is not useful for numerical calculation. 
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All Dq matrices exhibited in this paper satisfy rank(DQ) = 1 for some positive 
integer q. Further, there is an invertible m x m matrix Q such that Q~ 1 DqQ is the 
matrix whose (0, 0) th entry is 1 and all of whose other entries are 0. Define 

D' = Q~ X D Q Q, D[ = Q- 1 D 1 Q; 

by D' z (0,0) is meant the upper left corner entry of D' z . Let x(0 9 ) denote the set of 
all finite binary words z with no subwords q and with rightmost digit 1. 
Moshe [3] proved that 

1 i 

A = ^ r Y In ID' (0,0)1 

2^1(2^-1)^2^) 1 A n 

where £(z) is the length of z. We will prove in section |0.9j that 
2 ( t n 2 2q+1 - (3 + q)2 q + A - 



<r= 1 + 2 ^— f A^-2Ak + // 



where 



k=— 77^ r Y g l,ffl In ID' (0,0)1, 

0(7+1 (Oq _ -\\ ^ oi(z) 1 zy ' 71 ' 

z l z L ) ze x (0") z 

^= 2 ^-n E 2^) (ln| ^ (0 ' 0)l)2 - 

z ^ z ^ ze x (o«) z 

Summation of the series, coupled with Wynn's ^-process for accelerating convergence, 
serves as our primary method for calculating a 2 . For now, we revisit number-theoretic 
functions in [1] and compute both "average dispersion parameters" L(2)/ln(2) and 
"typical dispersion parameters" <7 2 /ln(2). 

0.1. Binomials. Define f(ri) to be the number of odd coefficients in (l + x) n . Let 
N denote a uniform random integer between and n — 1. We have q — 1, D — (1), 
Di = (2), 

D ®D = (1), D X ®D 1 = {A), e L W = 5/2 



ll6IlC6 

ln(Var(/(iV))) ln(5/2) 
ln(n) " ln(2) 

as n -> 00. Also D' Q = (1), Di = (2), 



1.3219280948873623478703194.. 



= 0.3465735902799726547086160..., 
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K 



1 ~ fc + 1 

- > — ; — In 

4^ n 2 k 

k=0 



1 00 1 

k=0 



ln(2) ~ fc(fc + l) 
4 ^ 

_ ln(2) 



2 k 



k=0 

2 00 1 2 
x - ft/ 

K=0 



21n(2), 



ln(2) 2 , 



cr 2 = 3A 2 -2Ak + /i = f--2 + -Jln(2) 
= 0.1201132534795503561667756... 



ln(2) ; 



hence 



Var(ln(/(JV))) ln(2) 



0.1732867951399863273543080... 



ln(n) 4 

as n — > 00, confirming a result of Kirschenhofer [4]. The typical parameter 0.173... 
is considerably smaller than the average parameter 1.321... because outlying values 
(which occur rarely) have been damped by the logarithm. 

0.2. Trinomials I. Define g(n) to be the number of odd coefficients in (1 + x + 
x 2 ) n . We have 5 = 1, 



D 



1 2 




D 1 



1 2 
1 



D ® D 



( 1 2 2 4\ 





V J 



hence 



ln(Var((/(JV))) ln(0 



Di®Di 



^ 1 2 2 4 ^ 

1 2 0* 

12 

V 1 / 



1.4924205743549514375202537... 



ln(n) ln(2) 

where £ = e^ 2 ** = 2.813... has minimal polynomial £ 3 — 2£ 2 — 3£ + 2. Also [5 J 



D' 



1 




D[ 



3 -4 
1 -2 



-r 



. 1 ~ 1 . /2 fc + 2 - ( 



l~Jfe + 1 /2 fc+2 -(-l) fc> 
k = - > — ; — In 



0.4299474333424527201146970. 



1™ 1, /2 fc + 2 
^ fc=o z v 



\k\ 2 
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and 



a 2 = 3A 2 - 2A k + n = 0.1211367118847285164803949... 
Var(ln(#(iV))) a 2 



ln(n) 



ln(2) 



0.1747633335056929866262498... 



as n — > oo. 



0.3. Quadrinomials. Define gz{ri) to be the number of odd coefficients in (1 + 
x + x 2 + x 3 ) n . This extends our earlier definitions / = g± and g = gi- We have q = 2, 





/ l 


2 






f o 







A) = 








;)• 




2 





i) 













v0 


1 





hence e L(2) = 5/2. Also 





f 1 








D'o = 













v0 


1 
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-4 


-6 





2 


1 


2 


-4 


-4 



hence pQ 



1 



A = 

2 + £(z 



E 



24 ^ 2'(*) 
- 4 z& x (00) z 



In |^(0, 0)| 



ln(2) 



Z4 z e x (oo) z 



^ 1^(0,0)1, „ = i £ -Lj 

Z4 ^Gx(OO) Z 



(ln|^(0,0)|) 5 



a 2 = 2|A 2 - 2A k + a* = 0.12011325..., 

and cr 2 /ln(2) = 0.17328679.... We conjecture that a 2 /ln(2) equals ln(2)/4 and prove 
this to be true in section |0.10j . 

0.4. Trinomials II. Define h 3 (n) to be the number of odd coefficients in (1 + x + 
x 3 ) n . This extends our earlier definition g = h 2 . We have q = 2, 



Dn 



( 1 


2 


1 


<M 




f 1 


1 


1 










1 


i 




1 








1 

















1 








V o 








o J 




V o 





1 


1 / 



hence 



ln(Var(/i 3 (iV))) ln(0 



ln(n) 



ln(2) 



1.5459492845008943975543991... 
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where £ = e L ^ = 2.919... has minimal polynomial 

16£ 10 - 40£ 9 - 36£ 8 + 22£ 7 + 76£ 6 + 7£ 5 - 19£ 4 - 19£ 3 + 2£ + 1. 

Also 



£>1 



^ 1 ^ 

10 

V o o o o / 

hence A = 0.45454538229305..., a 2 = 0.12497319... and 



(3 -6 -2 4 ^ 

1 10* 
1-3-2 2 

\01 00/ 



Var(ln(/i 3 (iV))) _ a 2 
Hn) * ln(2) 



= 0.18029820... 



as n — > oo. 



0.5. Quintinomials. Define ^(n) to be the number of odd coefficients in 
x + x 2 + x 3 + a; 4 )™. We have q — 2, 



/ 1 1 2 \ 



10 2 

V / 



/ 1 2 \ 

10 

10 2 

V 1 J 



hence 



ln(Var(<7 4 (A0)) ln(0 



1.6534827473445406557431504. 



ln(ra) ln(2) 
where £ = e L ^ = 3.145... has minimal polynomial 

4£ 10 - 8£ 9 - 21£ 8 + 14£ 7 - 28£ 6 + 126£ 5 + 65£ 4 + 68£ 3 + 48£ 2 - 56£ - 32. 

Also 



D> 



( 1 \ 

10 

V o o o o J 

hence A = 0.504253705692..., a 2 = 0.11406217... and 



/ 5 -10 -8 4 \ 

1 -1 -2 -2 

1-3-2 4 

V 1 -2 J 



VarQnMAQ)) 
ln(n) 



m(2) 



= 0.16455692... 



as n — > oo. 
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0.6. Trinomials III. Define h±{n) to be the number of odd coefficients in 
x + x 4 ) n . We have q — 2, 
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/ 



hence 



ln(Var(/i 4 (A0)) ln(0 



ln(n) 



ln(2) 



/ioiooioo\ 

10000000 
01011001 
01000000 
00100001 
00010010 
00001000 
\00000110/ 



1.5707744868006419128591802... 



where £ = e L ^ = 2.970... has minimal polynomial 



32£ 13 - 80£ 12 - 8£ n - 60£ 10 - 232£ 9 + 240£ 8 
+44£ 7 + 9£ 6 + 11£ 5 - 54£ 4 - 4£ 3 + 3£ 2 + £ + 2. 



Also 
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/ 



hence A = 0.45759385431410..., a 2 = 0.13055386... and 

Var(ln(/i 4 (JV))) 



a 



ln(n) 



ln(2) 



0.18834940... 



as n 



oo. 
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0.7. Sextinomials. Define g^(n) to be the number of odd coefficients in (f + x + 
. . . + x 4 + x 5 ) n . We have q = 3, 



Dr. 



hence 



/ 1 1 2 2 \ 



f f f 





V f J 



ln(Vax(0 5 (JV))) ln(0 



D l 



/000000\ 
2 2 

f f 2 2 
f 

\ / 



f.69037507596394449f 5537652. 



where £ = e 



L(2) 



ln(ra) ln(2) 
- 3.227... has minimal polynomial 

f 28f 11 - 640£ 10 + 41 6£ 9 + 1008£ 8 + 416£ 7 - 28£ 6 
-31 12£ 5 - 2572£ 4 + 346£ 3 + 1887£ 2 + 511f + 144. 



Also 



D' = 



( 1 














o \ 
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o J 



D[ 



( 6 


-8 
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-10 
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-6 \ 
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-4 


-4 
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-3 


V o 
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-2 


1 ) 



hence 



1 



1 



A = 7^ £ ^777 In |^(0, 0)| = 0.5344481528..., 

iiZ 2€x(000) 



1 



1 



and 



h= y 3+ ,y in i^(o, o)i, // — y — 

112 ^oo) 2^ A " 112^ 00) 2^) 

a 2 = ±f\ 2 -2\k + (i = 0.0965... 

Var(ln((? 5 (iV))) a 2 



(ln|^(0,0)|) 2 , 



ln(ra) 



ln(2) 



0.1392... 



as n — > oo. 



Typical Dispersion and Generalized Lyapunov Exponents 



8 



0.8. Septinomials. Define g&{n) to be the number of odd coefficients in (1 + x + 
• • • + x 5 + x 6 ) n . We have q = 3, 



hence 



/ 1 1 2 \ 



1 2 

2 10 10 



\0 0/ 



ln(Var(<? 6 (A0)) HO 



where £ = e 



L(2) 



Infra) ln(2) 
- 3.307... has minimal polynomial 



/ 2 1 \ 

1 

1 2 

2 1 

1 

V 1 J 



1.7258729504941114967801068. 



+ 4072£ 15 + 302£ 14 



8£ 21 - 4£ 20 - 18£ 19 - 335£ 18 + 34£ 17 + 474£ 16 
-3119£ 13 - 16848£ 12 - 1056f n + 7321£ 10 + 29681£ 9 + 910£ 8 
-6690£ 7 - 22628£ 6 - 152£ 5 + 1936£ 4 + 6112£ 3 - 128£ - 512. 



Also 



D> 



( 1 \ 



1 

1 



Voooooo; 



f 7 


-36 


-28 


-24 
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-4 \ 
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-12 
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V i 


-6 


-6 


-4 


1 


o J 



hence A = 0.53765282..., a 2 = 0.1082... and 

VarQnG^AQ)) 
ln(ra) 



ln(2) 



0.1561... 



as n 



oo. 



0.9. Proof of Formula for a 2 . Fix a nonnegative integer k. In the definitions 
of A and a 2 , we assumed that each binary word z of length k occurs with probability 
2~ k . This assumption is not necessary: let Pk(z) denote the (non-uniform) probability 
associated with z. Let 



Cfc(t) 



2 Pk( z ) 

£(z)=k 



\D, 



t real 
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then the radius of convergence of 



V(s,t) = J2 s e(z) Pe(z)(z) \\D Z ||* = E c k(t)s k , s complex 

z k=0 



IS 



lim c k 1/k = lim (e(||^||*)) 



-1/k 



-Lit) 



Hence Q(s,t) converges absolutely if \s\ < e~~ L ^ and diverges if \s\ > e~ L W. For 
convenience, we will write p(z) instead of pu z \ (z) from now on. 

We also postulated that there exists a positive integer q such that rank(DQ) = 1. 
This postulate can be weakened to the following: there exists a binary word z* for 
which rank(_D^») = 1. Any product of matrices involving D z * is also rank 1 since, if 
D z * = a (3 T , then 

D x D z *D y = (D x a)(p T D y ). 

Let x( z *) denote the set of all finite binary words w such that z* appears in the word 
w z* only at the end. Hence, given an arbitrary binary word z, it follows that either 



(0) 



z = w 

for some G x( z *) or 

z = (wWz*)(wWz*)(wV-Vz*) ■ ■ ■ (w^z*)(w^z*)w^ 

uniquely, where j > and £ x( z *) f° r an 111 the latter case, clearly 

\\D Z \\ = \(3 T D wW a\ ■ \(3 T D w( 2)a\ ■ ■ ■ \f3 T D w Q-r,a\ ■ \p T D w(j) a\ ■ (D w(oo )a)((3 T D w(0 )) 

Consider the words fixed and the index j as increasing; the final factor 

is thus immaterial. It follows that the convergence behavior of Q(s,t) is the same as 

Y J s^ {1) ^ P {w {l) z*W T D wW a\ t + 

Y, E s' (u ' (2)z * w(1)z * ) p(^ (2) ^^ (1) ^)l/3 T ^( 2 )«l i l/3 T J D w (i)«i t + 

u?( 2 ) tin 1 ) 

E E E ^ (3) ^ (2V ™^ + 

w (3) w (2) w (l) 

which equals F(s,t)/(1 — F(s,t)), where 



F(s,t)= Y s^ w ^p(wz*)\(3 T D v 

wex(z*) 



■a\\ 
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Therefore the smallest zero in modulus of 1 — F(s,t) determines the radius of con- 
vergence of Q(s,t) and the generalized Lyapunov exponent satisfies 

L(t) = -ln(s(t)) where F{s{t),t) = l. 

From the formula for L(t) in terms of s(t), we deduce that 

5(0) = 1 

since L(0) = 0, thus F(1,Q) = 1. By the Implicit Function Theorem, 



s'(0) 



dF(s,t)/dt 



dF(s,t)/ds 



Ft(l,0) 
: 1_ ^-(1,0) 



thus 

By the Chain Rule, 
s"(0) 



A = L'(0) 



d F t (s,t) 



s'(0) _ F t (l,0) 



s(0) F S (1,0) 



dtF s (s,t) 
d F t (s,t) 



s=l 
t=0 



s 

t=0 



dsF,(s,t) 
I F*t_ F ss F t 
\1\ 



x ■ m 



d F t (s,t) 



dtF s (s,t) 



s=l 
t=0 



^ 2 



thus 



a 2 = 



^(o) = - a(0) % ^ =^(0) a -^(0) 



It is easy to show that 



F t = E p(wz*)\n\p T D w a\, 

F st = E t(wz*)p(wz*)\n\f3 T D w a\, 

wex(z*) 

Ft= E ^(^^(inl^^al) 2 . 

wex(z*) 
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In our scenario, z* = 9 . The derivatives F s and F ss can be written in closed- form 
because 

X (0 q ) = {0} U xo U (xo x X o) U (xo x xo x xo) U • • • 

where Xo — {1, 01, 2 1, . . . , O 9-1 !} and the Cartesian product Xo x Xo is to be inter- 
preted as concatenation: 

{11,101,10 2 1,...,10 9-1 1, 
Oil, 0101, OIO^^.^OIO 9 - 1 !, 

Xo x xo = ° 211 ' ° 2l01 > ° 2l ° 21 ' • • • > O 2 ^" 1 !, 

9_1 11, 9-1 101, 9-1 10 2 l, . . . , 9_1 10 9_1 1} 

hence 



s \« - / 9 /s 



It follows that 



and therefore 



2/ V jfe! ^2/ / ^2/ V 2/ V V2 



F s (l,0) = 2(2 9 -1), 
F ss (l,0) = 4(2 29+1 -(3 + g)2 9 +l) 



c 2 = ^iA 2 + A 2 -2^A + ^ 
F F F 

A s A s A s 

( 2 29+1 - (3 + g)2 9 + 1 \ - F st F tt 

= ( 1+2 krr — j A - 2A if + l?' 

We conclude the proof by setting k = F st /F s and /x = F tt /F s . 

0.10. Proof of Conjecture. Before giving the proof, let us extend the formula 
for F(s,t) to the case of three nonnegative square matrices E , Ei, E 2 satisfying the 
following: there exist two ternary words z{, z\ for which rank(_E z *) = rank(i^) = 1. 
Write E z * = a,\ (3f, E z * = a 2 0£ and assume that z* is not a subword of z* for % ^ j. 
Let xi z h Z V) denote the set of all finite ternary words w such that z* appears in the 
word w z* only at the end, 1 < % < 2. A unique factorization property holds for all 
ternary words z as before. The function F(s, t) becomes a 2 x 2 matrix with entries 

= E s^^piwzDlPjE^. 



Typical Dispersion and Generalized Lyapunov Exponents 12 



The smallest zero in modulus of det(7 — F(s, £)) gives rise to the generalized Lyapunov 
exponent L(t) = — \n(s(t)) analogous to before. 

Let us return to binary words z. We assume without loss of generality that zn — 1 
(for this argument may be repeated with Os and Is interchanged). Every word thus 
looks like 

2 = 1 0*° 1 jl 10 j2 ... 1 o jn - 2 1 0**- 1 
where each j« > 0. We introduce a rewording of D z : 

D j0 D h D h ■ ■ ■!),., J) h , = (DM) (DM) (DM) • ■ • (dMo~ 2 ) (Aitf- 1 ) 

and note that the weight associated with D jz is (s/2) ji+1 . Also 



D = D 1 




D 1 = DiD Q 





Dj = DM 



for all j > 2. Hence the rewording actually consists of only three matrices D , Di, 
D 2 with weights s/2, s 2 /4 and 



oo / x m +l 
m=2 



S 3 A 



l-s/2 



Further, the initial row of each matrix is zero, thus we may consider only the lower- 
right 2x2 submatrix: 





1 2 

E 1 

4 4 ' 




(^(l 2)=a^, 



4 
1 



M, 



( J)(4 A)=a 2 pl 



Let ri(s) = s/2, r 2 (s) = s 2 /4 denote the weights corresponding to cti0[, a 2 P 2 and 
q(s) = (s 3 /8)/(l — s/2) denote the weight corresponding to M. Setting z x = and 
z* 2 = 2, we obtain 

oo 

^'(s,t) = 5:r,(s)g(s) fc |/5jM fc a 

fc=0 
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As calculated in pQ, 



2 if i = 1, j = 1 

2 2 if i = 1, j=2 

2 2k if i = 2, j = 1 

2 2(fc+i) if i = 2, j = 2 



hence 



*i,iM) = E 



00 s (s r 



fc=0 
00 



2 V 4 



2*C 

1 - c 2 
1-C 2 



00 3/ 



fc=0 

00 



s/2 V 4 



- 2 



2fct 



c 3 



^ = £1^72(4 



22(Jfc+l)i 



;i-0(i-2 2 *c 2 ) 

2 2 *C 3 

: (1_ C )(1_ 2 2* C 2) 



where ( = s/2. The smallest zero of 

(l-(2* + l)0(l-(l-2* + 2 2i K 2 ) 



det(F(s,t) -/) 



is ((t) = 1/(2* + 1), hence 



;i-0(i-c 2 )(i-2 2 *c 2 ) 



L(f) = -ln(2C(t))=ln((2« + l)/2). 

This expression for L(t) is trivially the same as the generalized Lyapunov exponent 
for binomials, thus A = ln(2)/2 and a 2 = ln(2) 2 /4. 

0.11. Digital Sums. Let #(n) denote the number of Is in the binary expansion 
of n. We know that #(n) = ln(/(n))/ln(2) and #(3n) = ln(p3(n))/ln(2), therefore 
#(iV) and #(3iV) are identically distributed in the sense that 



#(A0- 



ln(n) 
21n(2) 



\ 



V 



ln(n) 
41n(2) 



< X 



m I exp (4) * *- p 



#(3iV) - 



ln(n) 
21n(2) 



V 



ln(n) 
41n(2) 



< X 



/ 



as n — > 00. Moreover, #(aN + b) are identically distributed, < b < a. This answers 
a question raised in [1] and we refer interested readers to details in [6j Ej. 
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0.12. Closing Words. It is natural to seek a Central Limit Theorem for functions 
examined in this paper, for example, #(n) or f(n). This seems to be an open problem, 
but we indicate a possible direction for solution. Trollope [8] &; Delange [9] proved 
that 

nto 21,i(2) ' \te(2)) 

exactly, where <&(x) is a certain continuous nowhere-differentiable function of period 
1, 

-0.2075187496394219092731305... = - 1 = inf < sup$(x) = 0, 

2 ln(2) x ' x 



and the Fourier coefficients of are all known. The mean value of is [9], [12] 
i 

/ <5>(x)dx = — — ^ (ln(27r) - 1) - - = -0.1455994557083223046583226.... 
J V ' 21n(2) v v ; ; 4 

In principle, we can numerically compute any percentile of by approximating 
the Lebesgue measure of all x G [0, 1] satisfying $(a;) < r for some threshold r. This 
"inversion" must be done carefully, however, because of the fractal nature of <&(x). 
A plot of the limiting density function would be good to see! 
Stein [10J & Larcher [H] likewise proved that 



ln(2) 1 j ^ KM*)). 



where $?(x) is a continuous nowhere- mo notonic (but almost-every where differentiable) 
function of period 1, 

0.8125565590160063876948821... = inf < sup^(x) = 1, 

X X 

and the Fourier coefficients of $?(x) are all known. No closed-form expression is 
evident for the lower bound [13l [T41. The mean value of ^f(x) is 



^(x)dx = 0.8636049963990796049605033. 



and this too is unknown p^2l [15] . Again, a plot of the limiting density corresponding 
to ln(^(a;)) would be welcome progress. 

The function g(n) deserves more attention: Fourier expansions for the analogs of 
both $(a;) and ty(x) are desired (if these exist). We have not mentioned thus far the 
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functions u(n) associated with Pascal's rhombus or v(n) associated with "Fibonacci's 
rhombus" pQ. For Stern's sequence, in which v (n) is the number of odd coefficients 
in 

p n (x) = Xp n -i(x) +p n -2(x), px(x) = X, p Q (x) = 1 

we have 

A = 0.396212564297744..., 

a 2 = 0.022172945128737... 

and the method for calculating a 2 will be published later [16]. By contrast, u(n) is 
the number of odd coefficients in 

p n (x) = (1 + X + X 2 )p n _i(x) + X 2 p n - 2 {x), Pi (x) = 1 + X + X 2 , p (x) = 1. 

and A = 0.57331379313..., but no precise estimate of a 2 has yet been found. 
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